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Abstract: Countries with high lactase nonpersistence (LNP) or low lactase persistence (LP) 
populations have lower rates of some "western" diseases, mimicking the effects of sunshine 
and latitude. Inflammatory bowel disease (IBD), ie, Crohn's disease and ulcerative colitis, is 
putatively also influenced by sunshine. Recent availability of worldwide IBD rates and lactase 
distributions allows more extensive comparisons. The aim of this study was to evaluate the extent 
to which modern day lactase distributions interact with latitude, sunshine exposure, and IBD 
rates. National IBD rates, national distributions of LP/LNP, and population-weighted average 
national annual ultraviolet B exposure were obtained, estimated, or calculated from the literature. 
Negative binomial analysis was used to assess the relationship between the three parameters 
and IBD rates. Analyses for 55 countries were grouped in three geographic domains, ie, global, 
Europe, and non-Europe. In Europe, both latitude and ultraviolet B exposure correlate well with 
LP/LNP and IBD. In non-Europe, latitude and ultraviolet B exposure correlate weakly with LP/ 
LNP, but the latter retains a more robust correlation with IBD. In univariate analysis, latitude, 
ultraviolet B exposure, and LP/LNP all had significant relationships with IBD. Multivariate 
analysis showed that lactase distributions provided the best model of fit for IBD. The model 
of IBD reveals the evolutionary effects of the human lactase divide, and suggests that latitude, 
ultraviolet B exposure, and LP/LNP mimic each other because LP/LNP follows latitudinal direc- 
tions toward the equator. However, on a large scale, lactase patterns also follow lateral polarity. 
The effects of LP/LNP in disease are likely to involve complex interactions. 
Keywords: lactase, latitude, ultraviolet B exposure, evolution, inflammatory bowel disease 

Introduction 

The ability of adults to digest lactose in milk divides humanity into two phenotypes. Those 
able to digest this disaccharide in adulthood are considered lactase-persistent (LP), while 
those who lose this ability are called lactase-nonpersistent (LNP). The enzyme lactase 
phlorizin hydrolase, residing on the brush border of the proximal intestine, is geneti- 
cally determined. 1 Several polymorphisms have been identified which control lactase 
transcription in a cis position distally on chromosome 2. 1-5 There are distinct geographic 
patterns of the LP/LNP phenotype, with global population distributions of roughly one 
third LP and two thirds LNP. Modern day geographic distributions were determined 
by evolution and migrations 7,000-10,000 years ago (approximate time of appear- 
ance of LP phenotype). Independent emergence of the LP phenotype in Africa and the 
Middle East and further migrations after the discovery of the New World contributed 
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to modern day LP/LNP patterns. 6-8 The hypotheses regarding 
emergence of LP status concern latitude, sunshine, and ultra- 
violet B exposure (the calcium assimilation hypothesis) and 
simultaneous evolution of genes related to pastoralism and 
animal husbandry (the gene-culture coevolution hypothesis). 
In Africa and the Middle East, pastoralism is dominant, while 
both are possible in Europe. 1,9 ' 10 

Within the last six decades, studies have compared national 
lactase distributions with certain national disease rates. 1 1-13 In a 
previous publication, five "western" cancers and one "east- 
ern" cancer were compared according to distribution of LP/ 
LNP, and it was reported that cancers of the colon, ovary, 
prostate, breast, and lung were less common in populations 
with increasing LNP frequency, but stomach cancer was 
proportionately increased. These conclusions were based 
on the limited available data. 13 Crohn's disease (CD) and 
ulcerative colitis (UC) (both are types of inflammatory 
bowel disease [IBD]), common in industrialized nations, 
were found to share diminished rates with increased LNP 
status. 1213 In this same time period, diseases in industrialized 
regions have been observed to follow distinct geographic 
patterns, with a north-to-south gradient. These include 
cancers, 1415 hematological malignancies, 1617 multiple 
sclerosis, 18 CD, UC, 19,20-22 and other diseases. The puta- 
tive rate modifier in the north-to-south gradient paradigm 
may be sunshine, ultraviolet B exposure, and subsequently 
vitamin D synthesis. 1415 

Since the rate-reducing effects of lower latitude, higher 
ultraviolet B exposure, and higher frequency of LNP 
populations with some diseases are similar, we hypothesized 
that the combined effects of ultraviolet B exposure and 
LP/LNP distribution may have determined recent geographic 
disease patterns. 

In this paradigm, not only is there a north-to-south 
direction but also some lateral changes in regions. Recent 
availability of additional data on both IBD rates 23 and LNP 
frequencies 24 allow a re-evaluation of these diseases with 
respect to latitude, ultraviolet B exposure, and distribution 
of LP/LNP. In this context, IBD offers a model for evalua- 
tion of the hypothesis that LP/LNP serves an evolutionary 
background for a number of modern day diseases. 5 ' 25 

Materials and methods 

Literature search 

There are limited data available for matching regional IBD 
rates with regional lactase phenotypes. Therefore, national 
data were sought or estimated from the literature. Most data on 
disease rates for IBD were based on one recent publication. 23 



Some rates were obtained from other reviews. 26-29 PubMed 
and Google Scholar were also consulted to determine further 
incidence rates. For IBD, the MeSH words used were: 
"incidence, prevalence of Crohn's disease and/or ulcerative 
colitis around the world" or "international", "national rates of 
Crohn's disease", ulcerative colitis", "change in epidemiology 
or epidemiology" or "inflammatory bowel disease or Crohn's 
disease or ulcerative colitis". 

As a control against using average values for lati- 
tude and national IBD rates, CD rates for individual 
European cities as published by Molodecky et al were also 
correlated with their latitudes. 23 Latitudes for individual cit- 
ies were obtained from the Internet and are presumed to be 
more focused. 

Similarly, a recent review was used for lactase distribution, 24 
but older reviews were also searched. 1 1,3 °' 31 For further lactase 
distribution, the terms "lactase persistence or non-persis- 
tence", "lactose maldigestion or intolerance", or "genetics of 
lactase" were used. In each topic, individual references listed 
in relevant papers were also evaluated. Several authors were 
contacted directly for national estimates either on IBD rates 
(based on published reports 32 " 34 ) or LP proportions. 35 Finally, 
for lactase distributions, two online distribution maps were 
used as a rough estimate of LP/LNP proportions 24 ( http:// 
www.britannica.com/EBchecked/media/157598/Global- 
distribution-of-lactose-intolerance-in-humans ). 

Selection criteria 

As stated above, national rates were preferred, failing 
which regional rates that could be used to estimate national 
rates were sought. Several publications were also used in 
which rates could be deduced by the methods described 
above. These include intestinal biopsies 36 and genetic 
studies for lactase. 37 - 38 These were used to estimate LP/ 
LNP rates. 

Handling of data 

Nationwide rates were included for the time period described 
in the Results section. Because some of these were unavail- 
able as single values, available regional data for IBD or 
ethnic/racial distributions within countries for lactase propor- 
tions were used to calculate national rates, using equations 1 
and 2. Regional data were matched with populations around 
the time the data were published. Actual numbers of patients 
were calculated and added for each available region. Summed 
populations were calculated as representative of a nation 
matched for the population during that period. The total 
numbers were then proportioned to incidence rates in 10 5 . 
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Estimation of national disease rates (D) based on regional 
data where X. is the number of patients with new disease in 
region or city "i", A. is population of region/city, and P is 
the total population 

n 

I** 

n 

where P=XA; and n is the number of cities and/or 

i=l 

regions. 

Estimation of national LP or LNP rates (L) based on 
ethnic population percentage, E., and fractional ethnic 
population, f. 

m 

L = £f i E 1 [2] 

i=l 

where m is the total number of ethnic groups. 

National yearly ultraviolet B (280-315 nm) exposures 
were deduced from the data of Lee-Taylor and Madronich 39 
and have been described previously. 13 Briefly, monthly 
surface level radiation based on a radiative transfer model 
driven by satellite-measured variables was used. Annual 
averages from the sum of monthly averages for the period 
1990-2000 were computed. To obtain a single representative 
value for each of the countries, population-weighted averages 
for ultraviolet B surface radiation were calculated for the 
locations of the largest population centers in each country. 
A single population-weighted latitude was calculated for 
each country using the same population weighting as used 
for calculation of the population-weighted ultraviolet B 
exposure. 

Calculation of national annual average ultraviolet B 
exposure 

I^UVBj 

UVB = ^ [3] 

P 

where P are the populations of the N population centers 

1 N 

considered and here P = X Pi- 

i=l 

UVB. is the annual ultraviolet B exposure at population 
center i. The number of population centers (N) included in the 
calculation of a national average varied from one for small 
countries to typically ten or more for the larger countries 
with many large population centers. 

For the results that are presented in terms of a single 
latitude for each country, a population- weighted latitude was 



calculated for the country using the same population weight- 
ing as applied for calculation of the population-weighted 
ultraviolet B exposure. So, 

N 

S p . LAT i 

LAT = ^L _ [4] 

Slat 

where LAT. is the latitude of population center i. Local and 
national populations were obtained from census populations 
and corresponded to the median of the time range for the 
published observation periods. 

Statistical analysis 

Relationships between latitude, ultraviolet B exposure, 
LP/LNP distribution, and rates of CD and UC were 
assessed using Spearman correlation coefficients because 
of the non-normal distributions for most of the variables 
(especially IBD rates). The relationship between the three 
global variables and the IBD rates were first assessed 
using Spearman correlation coefficients, and then by 
univariate and multivariable analyses. Since correlations 
between some explanatory variables are very strong (eg, 
>0.9 between ultraviolet B exposure and latitude), when 
conducting the multivariable analysis, highly correlated 
variables (correlations >0.5) were not entered into the 
model at the same time. 

Both the dependent variables, ie, CD rate and UC rate, 
were positive integers showing an overdispersed Poisson 
distribution (variance greater than the mean). Hence, lin- 
ear models were not appropriate, and negative binomial 
regressions were used in both univariate and multivariable 
analyses. In a few cases was a negative binomial model did 
not converge at a specified limit (ie, the relative Hessian 
convergence criterion was greater than 0.0001), a Pois- 
son model that converged normally and provided similar 
estimates was used instead. The exponentiation of the 
parameter estimate (EPE) of the negative binomial regres- 
sion model can be interpreted in a multiplicative manner, 
or as a disease rate ratio. An EPE > 1 indicates that a higher 
value for a continuous independent variable is associated 
with a greater disease rate, and an EPE of <1 indicates 
the opposite. For example, an EPE of 0.7 for ultraviolet B 
exposure indicates that the disease rate will decrease by a 
factor of 0.7 when there is a 1 0-unit increase in ultraviolet B 
exposure, if everything else is the same with regard to all 
other covariates. 
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Analyses were carried out for three geographic areas, ie, 
global, Europe, and non-Europe. Sensitivity analyses were 
carried out for three scenarios. In two cases, the propor- 
tions of LNP for three countries (Argentina, Panama, and 
Malta) were based on line map data only, and as such were 
deemed less reliable. In the case of Oman, data were also 
obtained from a more accurate projected map referenced 
in Itan et al. 24 Thus, these data were handled in two ways. 
In scenario 1, the proportion of LNP was reduced by 50% 
for three countries (Argentina, Panama, Malta), while in 
scenario 2, the four countries were removed whenever 
LNP/LP was involved in the analysis. The third case 
involves the rate of CD in Hungary. The original analysis 
used a rate of 2.2 per 10 5 population based on Bernstein and 
Shanahan, 26 while calculations from Molodecky et al 23 based 
on two regions gave an estimate of 4.6 per 10 5 population. 
Thus, in scenario 3, the rate of CD in Hungary was taken 
as 4.6 instead of 2.2 per 10 5 population from the original 
analysis. 

As a control for the calculated data, we also analyzed CD 
rates in specific cities in Europe (from Molodecky et al 23 ) 
against city latitudes using Pearson's correlations. All statisti- 
cal analyses were performed using SAS version 9.3 software 
(SAS Institute Inc., Cary, NC, USA). 

Results 

Seventeen papers were used for IBD rates and 26 for lactase 
rates. These included three authors who responded to contact, 
ie, two for IBD incidence rates (Dr CC Figueroa from Chile, 
Dr I Hilmi from Malaysia) and one for lactase distribution 
(Dr J Rocha from Portugal). LP/LNP rates for Argentina, 
Panama, and Malta were based on an online distribution 
map as described earlier, and Oman was also derived from 
a map reference. 24 

Incidence rates were restricted to study periods between 
1980 and 2008. Whenever possible, the more recent data 
were used. Values for CD (38/52 countries, 73%) and UC 
(38/48 countries, 88%) were based on data published by 
Molodecky et al. 23 For CD and UC, respectively, 15.4% and 
23% of national rates were calculated. The IBD rates based 
on responses from Chile 32 and Malaysia 33 34 are listed as 
estimates (4% of total for each disease). 

National lactase phenotype distributions (n=54) were 
derived from a variety of published reports as outlined, 16/54 
(30%) were derived from Itan et al, 24 calculated lactase rates 
were 2%, while national rates for LP/LNP frequencies were 
estimated in 13%. 35 Rates for IBD and lactase distribution 
are listed in Table 1 . 



In general, indirect tests described for lactase distributions 
included intestinal biopsy and predominantly lactose toler- 
ance and lactose breath hydrogen tests. All three have been 
validated against lactase genotype. 40,41 It is also assumed that 
national lactase rates have changed more slowly over time, 
especially in the Old World. 

Data from 55 countries were obtained, with a few miss- 
ing data values. Countries with IBD rates, latitude, annual 
ultraviolet B exposure (kJ/m 2 ), and percent LNP are listed 
in Table 1. In subsequent statistical analysis, the measure 
of ultraviolet B exposure was further divided by 100 (ie, 
ultraviolet B exposure of 7,000 is referred to as 70) in order 
to avoid small estimates and to facilitate comparisons. The 
year 2000 was the median for populations (based on median 
years of observation) regarding national disease rates. As 
such, the data cover a period of about three decades. 

Table 2 shows Spearman correlation coefficients for all 
variables in the three groups, ie, global, Europe, and non- 
Europe. Examples of global correlations for LNP and CD 
and for LNP and UC are shown in Figures 1 and 2. Summary 
results for r values are shown graphically in Figure 3A-C. 
In Figure 3C, it can be seen that the influence of latitude 
and ultraviolet B exposure on LNP in non-Europe decreases 
dramatically compared with Figure 3A and B. Modest lati- 
tudinal effects in non-Europe are evident and less affected 
by ultraviolet B exposure. In effect, Figure 3A and B show 
inverse mirror images. In Figure 3C, the modest to moder- 
ate influence of LNP on CD and UC is retained in all three 
geographic groups. 

The results of univariate analyses for rates of CD and UC 
against the three independent variables (latitude, ultraviolet 
B exposure, and LNP) within the three geographical group- 
ings are shown in Table 3. For both rates, the three variables 
showed a significant relationship in all three groupings. 
Table 4 shows the models that provided the best fits from 
the multivariable analyses. The best-fit model was selected 
using the Akaike information criterion (a measure of the 
relative goodness of fit of statistical models). For CD, LNP 
dominated statistically, both globally and in Europe, while 
latitude shared dominance with LNP in non-Europe. For 
UC, LNP dominated statistically, both globally and in non- 
Europe, while latitude fit the best model in Europe. The values 
in Tables 3 and 4 are the EPEs from the negative binomial 
analyses, which can be interpreted in a multiplicative man- 
ner, or as a disease rate ratio. For example, the CD rate ratio 
is 0.78 for LNP in the global data; this means that for every 
10-unit increase in the value of LNP (ie, a 10% increase in 
LNP proportion), the CD rate is decreasing by 0.78 times, or 
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Table I National population-weighted latitudes and average annual calculated ultraviolet B exposures 
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24 
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68 


Mexico 


21 


1 02 


1 3,692 


70 




2.1 3* 


102.6 


2002 


47 


I | 


Poland 


51 


20 


4,01 1 


37.5 


0.1 


1.8 


38.2 


2003 


29 


69 


Portugal 


40 


8 


7,123 


40X 


2.99 


3.6 


9.9 


1990 


23 


37 


Romania 


46 


25 


5,367 


5SX 


0.5 


0.97 


21.8 


2002 


23 


29 


Saudi Arabia 


25 


45 


1 3,042 


53 


1.66 




22.3 


2003 


48 


24 


Serbia 


44 


21 


5,426 


35 


1.84* 


1.31* 


7.7 


1996 


23 


1 1 


Spain 


40 


4 


7,156 


34 


5.5 


8 


40.28 


2000 


23 


24 


Uruguay 


33 


56 


8,810 


65 


0.74 


4.26 


3.3 


2003 


49 


70 



Notes: Lactase nonpersistence frequencies are derived from the indexed references. Incidence/ 1 0 s of Crohn's disease and ulcerative colitis are derived from the indexed 
references. The population data were obtained from on-line information and are matched as closely as possible to median year of disease incidence acquisition. *Countries 
denoted where national rates are calculated; X, countries denoted where national rates are estimated; ** http.7/www.britannica.com/EBchecked/media/ 1 57598/Global- 
distribution-of-lactose-intolerance-in-humans ; ***calculation from Molodecky et al 23 for Hungary was 4.6. 

Abbreviations: CD, Crohn's disease; IBD, inflammatory bowel disease; inc, incidence; LNP, lactase nonpersistence; UC, ulcerative colitis; UVB, ultraviolet B exposure. 
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Table 2 Spearman correlation coefficients of Crohn's disease 
(CD), ulcerative colitis (UC), ultraviolet B exposure (UVB), latitude 
and lactase non persistence (LNP) 





Ulcerative 
colitis 


1 KID 

LINK 


1 l\/D 

UVB 


Latitude 


Global 










CD 


0.75 


-0.73 


-0.53 


0.56 


UC 


1 


-0.59 


-0.38 


0.44 


LNP 




1 


0.74 


-0.76 


UVB 






1 


-0.98 


Europe 










CD 


0.68 


-0.59 


-0.38 


0.39 




1 
1 


-0.42 


-0,41 


U.4 1 


LNP 




1 


0.74 


-0.71 


UVB 






1 


-0.99 


Non Europe 










CD 


0.79 


-0.54 


-0.40 


0.49 


UC 


1 


-0.68 


-0.35* 


0.55 


LNP 




1 


0.09* 


-0.22* 


UVB 






1 


-0.86 



Note: *AII values were statistically significant except those that have an * mark (ie, 
P-value >0.05). 



is 22% less. Therefore, populations with 50% LNP are esti- 
mated to have 22% lower CD incidence rates than populations 
with 40% LNP. The opposite effect is true for latitude. The 
UC rate ratio is 2.01 in Europe for latitude, which means that 
for every 10-unit increase in latitude (eg, latitude 50 versus 



40), the incidence of UC increases 2.01 times, or by 101%. 
Figure 4A-C shows the recorded global distributions of CD, 
LP, and LNP, and the incidence of UC. Demarcations are 
color-divided by quintiles. 

Sensitivity analysis 

In scenarios 1 and 2 (concerning LNP distribution), only 
non-European countries were affected, so analysis within 
Europe did not change under these two scenarios. Globally, 
correlation of LNP with other variables under both scenarios 
barely changed (all changes are less than 0.03). However, in 
non-European countries, all changes were small except for 
the correlation between LNP and CD, which changed from 
-0.54 to -0.41 when changing the LNP distributions of the 
three countries. Univariate analyses showed almost no change 
in the disease rate ratio (changes were no more than 0.01) for 
either scenario when compared with the original analyses. 
Multiple variable analyses, which concerned only CD in 
the non-European theater, showed small changes whereby 
the disease rate for LNP changed from 0.78 to 0.80 and for 
latitude from 1.58 to 1.74. 

Under scenario 3, whereby rate of CD in Hungary changed 
from 2.2 to 4.6 per 10 s population, only small changes were 
observed. The correlations between CD with all other vari- 
ables showed a maximum change of 0.03, while the disease 
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Figure I Graphic presentation of the global relationship between CD incidence rates expressed as log CD rate and distributions of LNP as national percentages. 
Abbreviations: CD, Crohn's disease; LNP, lactase nonpersistence. 
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Figure 2 Graphic presentation of the global relationship between UC incidence rates expressed as log UC rate and distributions of LNP as national percentages. 
Abbreviations: LNP, lactase nonpersistence; UC, ulcerative colitis. 



rate ratios showed a maximum change of 0.04. It is expected 
that changing these results affected the global theater much 
less than the European or non-European theaters because 
there were fewer countries for subterritorial analyses. 

Finally, to control for outcome in more focused areas, 
we also did a correlation between CD rates and latitudes in 
41 European cities (r=0.55, P=0.002, Table 5). These rate/ 
latitude relationships mimic those found based on extrapo- 
lated national correlations in Europe. Lactase distributions 
were not available for individual cities. The relationship with 
ultraviolet B exposure would be expected to be similar given 
that correlations between ultraviolet B exposure and latitude 
are very high in Europe. 

Discussion 

While the etiology of IBD remains largely unknown, multiple 
environmental factors have been implicated. 26,71 However, none 
have been confirmed as causative. Observations of the pro- 
gression of IBD indicate that UC generally precedes a rise 
in CD by about two decades. In western countries in general, 
UC incidence rates are stable or declining, while those of 
CD are increasing or starting to level off. 71 In developing 
countries and countries adopting the western lifestyle, rates 
of UC may be rising. Among the variables deemed to be 
important and relevant, has been the generally consistent 
finding that incidence and prevalence rates for IBD decrease 



in a north-to-south direction. The possible role of LP/LNP 
distribution among variables having an impact has received 
only occasional mention. 12 

In the current ecological evaluation, we note modest 
to moderate correlations between both forms of IBD and 
latitude and ultraviolet B exposure. The subanalysis of data 
from European cities with more focused geography and indi- 
vidually more homogeneous populations also supports the 
north-to-south gradient in Europe. CD and UC are moderately 
correlated but there are some differences observed between 
CD and UC relations. There is a slightly better correlation 
between CD and either latitude or ultraviolet B exposure in 
the global domain, but a slightly better correlation between 
UC with either latitude or ultraviolet B exposure in Europe. 
There are somewhat divergent correlations with latitude or 
ultraviolet B exposure in both forms of IBD in non-Europe, 
such that the correlation with UC is now not statistically 
significant. 

We also note that global correlations between 
latitude or ultraviolet B exposure and LNP (LP/LNP 
distribution) drops to very weak levels in non-Europe, 
suggesting that the global pattern is largely due to that 
found in Europe. The course of UC likely changed in 
three decades so that rates may have leveled off, while 
that of CD is increasing. Further, they reflect the roughly 
two-decade difference in behavior pattern between UC 
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Figure 3 Bar graph distributions of correlations of LNR CD, or UC with latitude (A), 
ultraviolet B exposure (B), and relationship of LNP with CD or UC (C) is shown in three 
geographic domains of global, Europe, or non-Europe. 

Notes: (A) and (B) are"mirror images" of the effects of latitude and ultraviolet B exposure 
on inflammatory bowel disease and lactase distributions. In non-Europe, correlations of 
LNP with either latitude or ultraviolet B exposure drop to weak correlations which 
are not statistically significant and the correlation of UC with ultraviolet B exposure 
becomes nonsignificant. (C) demonstrates that lactase distribution still correlates well 
with both forms of inflammatory bowel disease.This figure demonstrates an independent 
mechanistic effect of LP/LNP on inflammatory bowel disease which is different from the 
effects of latitude/ultraviolet B exposure. *Nonsignificant correlations (P>0.05). 
Abbreviations: CD, Crohn's disease; LNP, lactase nonpersistence; UC, ulcerative 
colitis; G, global; E, Europe; N, non-Europe; LP, lactase persistence. 

and CD. While migration and intermarriage change the 
phenotypes of lactose digestion (LP dominant), these 
are slower than changes in IBD rates. 23 As such, LP/ 
LNP distribution is likely to be more stable than disease 
distribution. 



The impact of LP/LNP on IBD remains robust in all 
domains including non-Europe and the multivariable analysis 
supports the primacy of LP/LNP. We interpret this as evidence 
that the effect of LNP (LP/LNP proportions) is independent 
from that of latitude and ultraviolet B exposure on IBD. 

Studies showing the north-to-south 
gradient effect on IBD rates 

Studies from the latter half of the last century generally 
observed that IBD rates diminish toward the equator, with 
rising rates in Australia and New Zealand. 71 This pattern 
of latitudinal change has also been reported from Europe 
as a whole 21,27 and within individual countries. The same 
gradient has been observed in northern France, 72 Scotland, 73 
and the USA. 74 75 In a study of the Nurses Health cohorts I 
and II, Khalili et al arbitrarily divided the USA into three 
latitudinal regions, and over the period of the study found 
that the highest latitudinal region was associated with the 
highest CD rates. 75 

No north-to-south gradient is found in Canada, largely 
because there are few inhabitants in the north. However, the 
west coast of the country has the lowest rates, while the east 
coast has the highest. 76 The north-to-south gradient for CD 
appears to be more robust than that for UC. For example, in 
a study by Nerich et al, no north-to-south gradient for UC 
was noted, but was noted for CD. 72 However, in Finland, no 
gradient for CD was reported, but a north-to-south gradient 
was reported for UC. Finland has the highest rates of UC 
reported to date. 77 To our knowledge, there is no detailed 
analysis from Australia and/or New Zealand, but higher 
rates would be expected in New Zealand, with a south-to- 
north gradient in Australia. Part of the observed gradient 
discrepancies between CD and UC may relate to the fact 
that UC precedes CD by about 20 years. 78 

The current report incorporating three decades of disease 
trends supports the published findings. In addition, a very 
strong correlation between latitude and ultraviolet B expo- 
sure was observed, which in turn was correlated with IBD. 
This finding provides evidence for this relationship between 
ultraviolet B exposure and IBD rates. 19 ~ 21,79 ' 80 

IBD, latitude, and the ultraviolet B/ 
vitamin D hypothesis 

The main effect of latitude on IBD and other such geo- 
graphically patterned diseases is thought to be related to 
sunshine. The effect of ultraviolet B exposure is to increase 
production of vitamin D in the skin. 81,82 Evidence has 
emerged in the last 1 5 years that vitamin D modifies the 
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Table 3 Disease rate ratios (95% confidence intervals) from Univariate Negative Binomial Analyses (for every 10-unit change in each 
variable) 



Disease 


Territory 


Latitude 


UVB 


LNP 


Crohn's 


Global (n SI) 


1.51 (1.24-1.84) 


0.85 (0.79-0.92) 


0.78 (0.72-0.83) 




Europe (n 27) 


1.49 (1.07-2.07) 


0.83 (0.70-0.98) 


0.79 (0.70-0.89) 




Non Europe (n 24) 


2.48 (1. 50-1. 10) 


0.75 (0.63-0.89) 


0.73 (0.67-0.79) 


Ulcerative 


Global (n 47) 


1.45 (1.23-1.72) 


0.87 (0.81-0.94) 


0.80 (0.74-0.86) 


colitis 


Europe (n 24) 


2.01 (1.36-2.97) 


0.75 (0.62-0.91) 


0.80 (0.67-0.94) 




Non Europe (n 23) 


1.71 (1.23-2.38) 


0.84 (0.74-0.96) 


0.78 (0.72-0.84) 



Abbreviations: UVB, ultraviolet B exposure; LNP, lactase nonpersistence. 



outcome in many diseases, including IBD, through modu- 
lation of immunological mechanisms with anticancer and 
autoimmunity effects. 19,81,82 Yet the lack of vitamin D is not 
the putative cause of IBD and other such diseases. While 
other variables associated with latitude, such as climate, 
have also been identified as potentially relevant, 83 elements 
of the western lifestyle are thought to be causative. 26 Indeed, 
current publications on IBD rates show a rise in countries 
with traditionally high ultraviolet B exposure. Further, 
differences in IBD rates within regional populations have 
been found in high latitude northern 84 and equatorial or 
southern hemispheres 85,86 with more or less similar ultra- 
violet B exposure. 

Relationship between latitude/ultraviolet B 
exposure and lactase distribution 

In Europe, lack of sunshine and the consequent reduction of 
vitamin D synthesis in the skin has been postulated to lead to 
strong selection for LP (the calcium assimilation hypothesis). 87 
As outlined in the introduction, a counter hypothesis supported 
by several groups is that LP status evolved in conjunction 
with pastoralism and herding (the gene-culture coevolution 
hypothesis). 1,8,10 Distribution of lactase phenotypes in the 
Old World may then have depended both on selection and 
migrations from central Europe to northern Europe, Russia, 
and India to the east. 8,88,89 Migrations to the New World and 
South Pacific regions generally were made up of people with 



Table 4 


Disease rate ratios 


(95% confidence 


intervals) from 


Multivariable Negative Binomial Analyses (for 


every 10-unit 


change in 


each variable) 






Disease 


Global 


Europe 


Non-Europe 


Crohn's 


n=5l LNP 0.78 


n=27 LNP 0.79 


n=24 LNP 0.78 




(0.72-0.83) 


(0.70-0.89) 


(0.71-0.85) 








Latitude 1.55 








(1. 10-2.17) 


Ulcerative 


n=47 LNP 0.80 


n=24 Latitude 2.0 1 


n=23 LNP 0.78 


colitis 


(0.74-0.86) 


(1.36-2.97) 


(0.72-0.84) 



Abbreviation: LNP, lactase nonpersistence. 



LNP, and occurred prior to the emergence of LP genetic domi- 
nance. These hypotheses help to explain the more frequent 
LNP status of indigenous populations in North and South 
America as well as in New Zealand. The southern continent 
of Australia and New Zealand was populated much later by 
LP people from the British Isles. 

As seen in Figure 3A-C, the apparent impact of latitude/ 
ultraviolet B exposure is markedly diminished in non-Europe. 
This effect could be due to the method of assigning single 
values for each variable to large countries like the USA, 
People's Republic of China, and Australia. However, this 
unusual pattern was also reported previously based on both 
indirect tests 90 and genetic frequencies of LP. 9 

Relationship between disease 
and lactase phenotype 

If we accept that the current study confirms previous reports 
of a relationship of latitude and ultraviolet B exposure with 
IBD and confirms previously reported models of the partial 
relationship of these variables with LP/LNP proportions, 
then we should consider the relationship between LP/LNP 
and IBD. The apparently similar relationship between all 
three test variables globally and in Europe is challenged by 
the minimal effect of latitude and ultraviolet B exposure on 
LNP and weaker correlations between ultraviolet B exposure 
and IBD. Yet LNP retains its impact on IBD and is in fact 
more robust. This dichotomy together with the hypothesized 
different global dispersions of lactase phenotypes suggests 
an independent LP/LNP mechanism in addition to that of 
latitude/ultraviolet B exposure. 

It is doubtful that LP distributions exert direct pathogenic 
effects. Rates of IBD have also risen in high LNP popula- 
tions (eg, Japan, People's Republic of China, and Korea). It 
is postulated that any modifier effect of an LP/LNP disease 
interaction would be more complex than that currently 
attributed to latitude/ultraviolet B exposure. Several pos- 
sible modifiers of IBD can be hypothesized and are related 
to digestion of lactose and its evolution. 
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Figure 4 Maps of the world showing distributions of incidence of quintiles of inflammatory bowel disease, ie, Crohn's disease, (A) ulcerative colitis (C) and percentage 
lactase persistence (B). 

Notes: The map of inflammatory bowel disease covers a time span of about three decades. Color patterns are from highest to lowest frequencies. Quintiles for Crohn's 
disease each represent countries with calculated estimate rates of incidence spanning multiples of 2.75/ 1 0 s cases (overall range 0-13.75), ulcerative colitis 3.48/ 1 0 s cases 
(overall range 1-24.8) and lactase persistence 18% of the population (overall range 4-92). Epidemiologically, ulcerative colitis precedes the emergence of Crohn's disease 
by about two decades, and in general ulcerative colitis rates are stabilizing in older western societies, but may be rising in regions which are adopting western/industrialized 
lifestyles. Relatively, LP/LNP distributions may change more slowly. This time difference in disease progression may account for some of the variability in the relationship 
of ulcerative colitis with LP/LNP. Latitudinal changes in inflammatory bowel disease have also been described in individual countries. Furthermore, it is noted that lactase 
distributions in (B) show in large scale both latitudinal reduction in LP toward the equator (LNP) as well as lateral changes from politically western toward eastern nations. 
Abbreviations: LP, lactase persistence; LNP, lactase nonpersistence. 
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Table 5 Comparison of latitudes of European cities with incidence 

City Latitude CD inc/10 5 Population 

per 10 s 



Notes: Populations are based on median from the year 2000. Adapted from 
Gastroenterology, 142, Molodecky NA, Soon IS, Rabi D, et al, Increasing incidence 
and prevalence of the inflammatory bowel diseases with time, based on a systematic 
review, 46-54, copyright (2012), with permission from Elsevier. 23 
Abbreviations: CD, Crohn's disease; inc, incidence. 

Possible modifiers of IBD 
by lactase distribution 

Dairy food consumption 

In general, possible modifier mechanisms for LP/LNP sta- 
tus could relate to dairy food consumption. In this instance, 
a harmful effect of dairy foods for some diseases, including 
IBD, may occur and LNP status (by virtue of symptomatic 
lactose intolerance and culture) would reduce the amount 



of dairy foods consumed. In the case of IBD, a modest 
increased odds ratio was found on epidemiological grounds 
and was significant only for UC. 13 In another situation, LNP 
may be protective despite regular dairy food consumption 
via the effects of undigested lactose on microbial flora. 
In such cases, lactose could act as a prebiotic, promoting 
protective bacteria like bifidobacteria. 91 This scenario may 
represent an ecological fallacy on comparison of outcomes 
in epidemiological and patient-based studies, as has been 
shown in the case of colorectal cancer. 13 92 While mecha- 
nisms pertaining to this issue are controversial, dairy food 
consumption in high LNP countries is much less common 
than in high LP countries. 13 Nevertheless, consumption of 
dairy foods has been implicated in promoting development 
of IBD, 93 but further studies are required because outcomes 
are inconsistent. For example, a British study suggested a 
possible protective effect of unpasteurized milk consump- 
tion in CD. 94 

Coevolutionary genetics 

Coevolution of other genes with LP status could predispose 
to IBD in modern environments. Emergence of LP status is 
hypothesized to have incurred an increased risk of intestinal 
infections and death as a result of drinking unpasteurized 
milk. A number of genes may have coevolved with LP. For 
example, HLA types, leading to immune signaling, 95 and the 
NOD2/CARD15 system are hypothesized to have evolved 
in response to drinking milk and the threat of infection. 96 
Interestingly, genetic mutations in the NOD2 system, which 
are associated mainly with terminal ileal CD, are largely 
confined to Caucasians. In vitro at least, normal function of 
non-mutated NOD2 requires adequate l,25(OH) 2 D 3 to be 
present in cell culture medium. However, mutations inNOD2 
cannot be overcome with increasing vitamin D. 97 This in vitro 
study by Wang et al potentially links vitamin D with IBD, and 
may provide a plausible explanation as to how low ultraviolet 
B exposure may promote IBD at high latitude. 97 

Cystic fibrosis, which in the homozygous state is 
a lethal disease in Europeans, is hypothesized to have 
incurred a selective advantage in limiting diarrhea in the 
heterozygous form. 98 However, reports of possible protec- 
tive effects of heterozygous mutations against IBD are 
conflicting. 99100 

Overall, there are only a few contradictory studies 
examining the impact of lactase phenotype on IBD. 101-104 
However, almost two dozen other immune-related diseases 
share genotypes with IBD, 105 possibly linking a common 
epidemiology. 
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Possible modification of IBD by influence 
of geography on socioeconomic factors 

Countries with higher LNP frequency populations tend 
to be economically disadvantaged. This is true largely in 
South America, Africa, and south Asia. Countries like the 
People's Republic of China, South Korea, and Japan have 
already or are adopting more western lifestyles. Most of 
these countries had lower rates of IBD until towards the end 
of the last century. It is therefore possible that reporting of 
IBD rates is influenced by lack of experience with IBD or in 
economically less favored countries by lack of doctors and 
inexperience with the disease. 

In the middle of the last century IBDs had higher mortality 
rates. UC initially had a six-fold mortality rate compared with 
CD, but in the two decades following, the mortality of CD 
increased. 78 If such mortality rates were encountered in coun- 
tries with previously low IBD rates, their reporting would likely 
have been more thorough. However, no increased mortality 
has been reported in low IBD-incidence countries. This could 
suggest that mortality of IBD has spontaneously receded or 
that with reduced mortality these diseases are less recognized 
by local medical communities. 

Initial observations suggested that IBD was more com- 
mon in higher socioeconomic groups. Indeed, the study 
by Nerich et al from northern France did show a negative 
association between CD and farming and households below 
the poverty level. 72 However, other reports have refuted this 
claim. 106 Along similar lines, deprivation scores were not 
related to juvenile onset of CD in Scotland. 73 Therefore, the 
relationship has not been settled. It is possible that failure to 
diagnose or failure to get medical assistance in low socio- 
economic countries nevertheless contribute to low reported 
incidence rates. The appropriate reporting of outcomes of 
diseases more common in higher LNP populated countries 
(eg, stomach cancer, nasopharyngeal cancer, and hepatocel- 
lular cancer) suggests that unfamiliarity with IBD might 
contribute somewhat to under-reporting. More research on 
the impact of socioeconomic parameters on IBD rates would 
need to be carried out. 

Other possible effects on IBD 

Independent of dairy foods and genetic predisposition, 
observations that infectious diseases are still rampant in areas 
inhabited by large LNP populations 107 raise the question of 
their role in infections modifying host immunity. Immune 
reactions to such agents have been postulated to protect 
against modern day allergic and autoimmune diseases, 
including IBD 108 



Limitations 

A limitation of this analysis is the lack of ability to match 
data for the five variables (ie, latitude, ultraviolet B exposure, 
LP/LNP distribution, CD, and UC) in specific regions. As 
a consequence, it was felt that the best match would be to 
extrapolate data to a national level. Another limitation of 
this method is that large countries were identified by unique 
values for each of the five variables. This technique could 
hide existing relationships, such as the reports of north-to- 
south gradients in countries like the USA 74,75 and France. 72 For 
example, the nonsignificant latitude/ultraviolet B exposure 
relationship with LP/LNP and UC in non-Europe may reflect 
this potential problem. A third limitation is the wide time 
period chosen to evaluate IBD rates. The primary reason for 
this is the variable time periods for IBD rates reported from 
different locales. 

However, on examining the five evaluated variables, it 
should be noted that latitudes and ultraviolet B exposure are 
stable over time. Single values used for national descrip- 
tions take into account measurements from different parts 
of countries and are population-weighted. These unique 
descriptors do have a relationship to "national average" yearly 
ultraviolet B exposure and latitude. In addition, ultraviolet B 
exposure is more independent of the direction of latitude 
north or south of the equator. 

Although national lactase distributions are even less 
well defined, they are potential estimates of percentages of 
populations in local regions and may represent the majority 
of the country. However, in general, lactase distributions 
are more stable in Asia, Africa, and Australia, and to some 
extent in Europe. The most frequent disparity within a 
region would occur in Africa where tribal differences exist, 
but for our purposes, no matching IBD data are reported in 
such countries (eg, South Africa). It is true that migrations 
in the last decade and a half might have changed the land- 
scape somewhat, but more North Americans than Asians 
are able to digest lactose. This large population split results 
in lateral polarity outside Europe. The rates of change in 
IBD are likely greater than changes in national lactase 
proportions. 23 As such, disease rates are the predominant 
altering variables. The effect of these changes in disease 
rate is to reduce correlative findings because of stabilizing 
UC rates and expanding CD rates in different countries. 
The different times of onset and change in incidence rates 
between the two forms of IBD could be separated by as 
much as 20 years. 78 This time differential could account 
for the somewhat different relationships of CD or UC with 
LP/LNP distributions. 
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Strengths 

Our original hypothesis was that the LP/LNP evolutionary 
divide millennia ago has continued to have an impact on the 
geographic distribution of modern diseases in the last six 
decades. As a result, a general pattern correlation may be 
sufficient without precise regional effects. Within-country 
relationships will require far more appropriate data to be 
available in multiple regions of the world. Our findings are 
based on data that were independently published from dif- 
ferent sources, free of bias. 

We also emphasize in the current paradigm that neither 
the concept of the disease rate-reducing effects of latitude/ 
ultraviolet B exposure or LNP population frequency nor the 
relationship between LNP and latitude/ultraviolet B exposure 
are new. The most important finding here is that for the first 
time we relate these three independent variables simultane- 
ously to rates of IBD. 

The precise effect of lactase status on IBD and 
other diseases is not obvious (as was perceived for 
ultraviolet B exposure and vitamin D), but it is difficult to 
ignore when evaluating different disease rates among differ- 
ent racial groups (eg, first nations in North America, African 
Americans in North America, and Maoris in New Zealand 
which are all predominantly comprised of LNP people). 

Conclusion 

This analysis relates modern lactase phenotype distribu- 
tions to patterns of geographic distributions for IBD in 
the last 60 years. In the process, the previously observed 
relationships of latitude and IBD are supported. The close 
correlation between ultraviolet B exposure and latitude is 
confirmed and lends support to the ultraviolet B exposure/ 
vitamin D hypothesis and its impact on IBD. Similarly, the 
previously hypothesized and observed relationships between 
geographic distribution of LP/LNP and latitude/ultraviolet B 
exposure are supported. The findings of this study support 
the conclusion that the evolution of LP together with pre and 
post divide migrations contributed to the observed geographic 
spread of modern diseases currently common in western/ 
industrialized societies. 

We propose that the relationship between latitudinal 
polarity (north-to-south, south-to-north) and LP/LNP distri- 
butions mimic each other's disease-modifying effects through 
global LP-LNP polarity on a large scale as well. This is 
most clear in Europe, but there is polarity in non-Europe as 
well, where LP-LNP polarity generally exists from north- 
to-south or south-to-north, but seemingly less robustly so 
than in Europe. Further, polarity exists more evidently in 



lateral directions, and the current analysis shows that LP/LNP 
distributions can have independent effects from ultraviolet B 
exposure. We do not use the term "longitudinal" because 
this requires reference points; however, in geographic terms, 
there is a west-to-east and east-to-west polarity implied. In 
Europe, the polarity appears to be restricted to north-to-south 
and west-to-east. 

IBD represents a model for other diseases associated with 
polar distributions. The implications of the effects of LP/LNP 
are likely complex and would require investigation, perhaps 
for each related disease and on multiple levels. However, 
coevolution of genes predisposing to modern disease with 
evolution of lactase is at least one such likely contributing 
variable. 

Acknowledgments 

We would like to thank Dr Ian Shrier, Department of 
Epidemiology, Jewish General Hospital, McGill University, 
for his insightful advice and help in setting up data. We would 
also like to thank Professor Brian Smith, Desautel Department 
of Management, McGill University, for reviewing and com- 
menting on the manuscript. Finally, we thank Ms Melissa 
Szilagyi for assistance with calculations of national disease 
and lactase distribution rates. 

Disclosure 

This work was supported by a two-year stipend awarded to 
AS by the Department of Medicine, Jewish General Hospital. 
AS receives nonspecific general research funds from the 
Department of Medicine, Jewish General Hospital. In the past 
2 years, AS has participated in an advisory board meeting 
for the Dairy Bureau of Canada. The other authors have no 
conflicts of interest to declare. 



References 

1. Swallow DM. Genetics of lactase persistence and lactose intolerance. 
Ann Hum Genet. 2003;37(4): 197-219. 

2. Kruse TA, Bolund L, Grzeschik KH, et al. The human lactase-phlorizin 
gene is located on chromosome 2. FEBS Lett. 1988;240(l-2): 123-126. 

3. EnattahNS , SahiT, Savilahti E, et al. Identification of a variant associated 
with adult-type hypolactasia. Nat Genet. 2002;30(2):233-237. 

4. Tishkoff SA, Reed FA, Ranciaro A, et al. Convergent adaptation of human 
lactase persistence in Africa and Europe. Nat Genet. 2007;39(1):31— 40. 

5. Jones BL, Swallow DM. The impact of cis-acting polymorphisms on the 
human phenotype. Hugo J. 201 1;5(1 — 4): 13-23. 

6. Forster R Matsumura S. Did early humans go north or south? Science. 
2005;308(5724):965-966. 

7. Burger J, Kirchner M, Bramanti B, Haak W, Thomas MG. Absence of 
the lactase-persistence-associated allele in early Neolithic Europeans. 
Proc Natl Acad Sci USA. 2007;104(10):3736-3741. 

8. Itan Y, Powell A, Beaumont MA, Burger J, Thomas MG. The origins 
of lactase persistence in Europe. PLoS One Comput Biol. 2009;5(8): 
el 000491. 



Clinical Epidemiology 2014:6 



submit your manuscript | www.dovepress.com 
Dovepress 



195 



Szilagyi et al 



Dovepress 



9. Gerbault P, Moret C, Currat M, Sanchez-Mazas A. Impact of selection 
and demography on the diffusion of lactase persistence. PLoS One. 
2009;4(7):e6369. 

10. Simoons FJ. Primary adult lactose intolerance and the milking habit: 
a problem in biologic and cultural interrelations. II. A culture historical 
hypothesis. Am J Dig Dis. 1 970; 1 5(8):695-7 1 0. 

1 1 . Cramer DW, Xu H. Lactase persistence, galactose metabolism and milk 
consumption as risk factors for ovarian cancer. In: Aurichio S, Semenza G, 
editors. Common Food Intolerances 2: Milk in Human Nutrition and 
Adult-Type Hypolactasia. Basel, Switzerland: Karger; 1993. 

12. Nanji AA, Denardi FG. Primary adult lactose intolerance protects 
against development of inflammatory bowel disease. Med Hypotheses. 
1986;19(l):l-6. 

13. Shrier I, Szilagyi A, Correa JA. Impact of lactose containing foods and 
the genetics of lactase on diseases: an analytical review of population 
data. Nittr Cancer. 2008;60(3):292-300. 

14. Grant WB. A meta-analysis of second cancers after a diagnosis of 
non-melanoma skin cancer: additional evidence that solar ultraviolet -B 
irradiance reduces the risk of internal cancers. J Steroid Biochem Mol 
Biol. 2007;103(3-5):668-674. 

15. Grant WB. How strong is the evidence that solar ultraviolet B and 
vitamin D reduce the risk of cancer? Dermatoendocrinol . 2009; 1(1): 
17-24. 

16. van Leeuwen MT, Turner JJ, Falster MO, et al. Latitude gradients 
for lymphoid neoplasm subtypes in Australia support an associa- 
tion with ultraviolet radiation exposure. Int J Cancer. 20 13; 133(4): 
944-951. 

17. Mohr SB, Garland CF, Gorham ED, Grant WB, Garland FC. Ultraviolet 
light B and incidence rates of leukemia worldwide. Am J Prev Med. 
2011;41(l):68-74. 

1 8. Simpson S Jr, Blizzard L, Otahal P, Van der Mei I, Taylor B. Latitude is 
significantly associated with the prevalence of multiple sclerosis: a meta- 
analysis. J Neurol Neurosurg Psychiatry. 201 1;82( 10): 1 132-1 141. 

19. Ananathakrishnan AN, Khalili H, Higuchi LM, et al. Higher predicted 
vitamin D status is associated with reduced risk of Crohn's disease. 
Gastroenterology. 2012;142(3):482^189. 

20. Peyrin-Biroulet L, Oussalah A, Bigard M-A. Crohn's disease: the hot 
hypothesis. Med Hypotheses. 2009;73(l):94-96. 

2 1 . Shivinanda S, Lennard- Jones J, Logan R, et al. Incidence of inflamma- 
tory bowel disease across Europe: is there a difference between north and 
south ? Results of the European Collaborative Study on Inflammatory 
Bowel Disease (EC-IBD). Gut. 1996;39(5):690-697. 

22. Juste RA. Crohn's disease and ruminant farming. Got lactase? Med 
Hypotheses. 2010;75(1):7-13. 

23. Molodecky NA, Soon IS, Rabi D, et al. Increasing incidence and 
prevalence of the inflammatory bowel diseases with time, based on a 
systematic review. Gastroenterology. 2012; 142( 1 ):46— 54. 

24. Itan Y, Jones BL, Ingram C JE, Swallow DM, Thomas MG. A worldwide 
correlation of lactase persistence phenotype and genotypes. BMCEvol 
Biol. 2010;10:36. 

25. Szilagyi A. Determinants of geographic patterns of diseases: interaction 
of lactose/lactase status and sunshine exposure. Med Hypotheses. 
2010;75(5):466^170. 

26. Bernstein CN, Shanahan F. Disorders of a modern lifestyle: reconciling 
the epidemiology of inflammatory bowel diseases. Gut. 2008;57(9): 
1185-1191. 

27. Economou M, Pappas G. New global map of Crohn's disease: genetic, 
environmental, and socioeconomic correlations. Inflamm Bowel Dis. 
2008;14(5):709-720. 

28. Thia KT, Loftus EV, Sandborn WJ, Yang S-K. An update on 
the epidemiology of inflammatory bowel disease in Asia. Am J 
Gastroenterol. 2008; 103( 12):3 167-3 1 82. 

29. Lakatos L, Lakatos PL. Is the incidence and prevalence of inflamma- 
tory bowel diseases increasing in Eastern Europe? Postgrad Med J. 
2006;82(967):332-337. 

30. Sahi T. Genetics and epidemiology of adult-type hypolactasia. Scand 
J Gastroenterol. 1994;29 Suppl 202:S7-S20. 



3 1 . Scrimshaw N, Murray EB. The acceptability of milk and milk products 
in populations with a high prevalence of lac tose intolerance. Am J Clin 
Nutr. 1988;48 Suppl 4:S1083-S1 159. 

32. Figueroa CC, Quera PR, Valenzuela EJ, Jensen BC. Enfermedades 
inflamatorias intestinales: Experiencia de dos centres chilenos. [Inflam- 
matory bowel disease: experience of two Chilean centers]. Rev Med 
CM. 2005;133(11):1295-1304. Spanish. 

33. Hilmi I, Singh R, Ganesananthan S, et al. Demography and clini- 
cal course of ulcerative colitis in a multiracial Asian population: 
a nationwide study from Malaysia. J Dig Dis. 2009;10(l):15-20. 

34. Chua KH, Hilmi I, Ng CC, et al. Identification of NOD2/CARD 1 5 
mutations in Malaysian patients with Crohn's disease. J Dig Dis. 
2009; 10(2): 124-130. 

35. Coelho M, Luiselli Bertorelle G, Lopes Al, Seixas S, Destro-Bisol G, 
Rocha J. Microsatellite variation and evolution of human lactase 
persistence. Hum Genet. 2005; 1 17(4):329-339. 

36. Blomme B, Gerlo E, Hauser B, Vandenplas Y. Disaccharidase activities 
in Belgian children: reference intervals and comparison with non-Belgian 
Caucasian children. Acta Pediatr. 20O3;92(7):806-810. 

37. Timpson NJ, Brennan P, Gaborieau V, et al. Can lactase persistence 
genotype be used to reassess the relationship between renal cell car- 
cinoma and milk drinking? Potentials and problems in the application 
of Mendelian randomization. Cancer Epidemiol Biomarkers Prev. 
2010;19(5):1341-1348. 

38. Weinberg RB. Apolipoprotein A-IV-2 allele: association of its 
worldwide distribution with adult persistence of lactase and specu- 
lation on its function and origin. Genet Epidemiol. 1999; 17(4): 
285-297. 

39. Lee-Taylor J, Madronich S. Climatology ofUV-A, UV-B, and Erythemal 
Radiation at the Earth 's Surface, 1979-2000. Technical Note TN- 
474+STR. Boulder, CO, USA: National Center for Atmospheric 
Research; 2007. 

40. Ennattah NS, Kuokkanen M, Forsblom C, et al. Correlation of intestinal 
disaccharidase activities with the C/T-13910 variant and age. World J 
Gastroenterol. 2007;13(25):3508-3512. 

41. Marton A, Xue X, Szilagyi A. Meta-analysis: the diagnostic accuracy 
of lactose breath hydrogen or lactose tolerance tests for predicting the 
North European lactase polymorphism C/T- 13910. Aliment Pharmacol 
Ther. 2012;35(4):429^140. 

42. Vahedi H, Merat S, Momtahen S, et al. Epidemiologic characteristics 
of 500 patients with inflammatory bowel disease in Iran studied from 
2004 through 2007. Arch Iran Med. 2009;12(5):454-460. 

43. Zouiten-Mekki L, Zaouali H, Boubaker J, et al. Cardl5/Nod2 in a 
Tunisian population with Crohn's disease. Dig Dis Sci. 2005;50(1): 
130-135. 

44. Petritsch W, Fuchs S, Berghold A, et al. Incidence of inflammatory 
bowel disease in the province of Styria, Austria, from 1997 to 2007: 
a population-based study. J Crohns Colitis. 20 1 3;7( 1 ):58— 69. 

45 . Jussila A, Virta LJ, Kautiainen H, et al. Increasing incidence of inflam- 
matory bowel diseases between 2000 and 2007: a nationwide register 
study in Finland. Inflamm Bowel Dis. 2012;18(3):555-561. 

46. Kiudelis G, Zvirbytiene A, Sventoraitye J, Kupecnikas L. Referral 
epidemiology of ulcerative colitis and Crohn's disease in Kaunas 
region Lithuania: one year prospective study. J Crohns Colitis. 2008; 
2 Suppl 1:S70. 

47. Yammamoto-Furusho JK. Clinical epidemiology of ulcerative colitis in 
Mexico: a single hospital-based study in a 20-year period ( 1 987-2006). 
J Clin Gastroenterol. 2009;43(3):221-224. 

48. Fadda MA, Peedikayil MC, Kagevi I, Kahtani KA, Ben AA, Al HI. 
Inflammatory bowel disease in Saudi Arabia: a hospital-based clinical 
study of 312 patients. Ann Saudi Med. 2012;32(3):276-282. 

49. Buenavida G, Casanias A, Vasquez C, et al. National Inflammatory 
Bowel Disease Registry. Acta Gastroenterol Latinoam. 201 1 ;4 1 (4): 
281-287. 

50. Wang YG, Yan YS, Xu JJ, Du RF, Kuhnau W, Flatz G. Prevalence of 
primary adult lactose malabsorption in three populations of northern 
China. Hum Genet. 1984;67(1): 103-106. 



submit your manuscript | www.dovepress.com 
Dovepress 



Clinical Epidemiology 2014:6 



Dovepress 



Lactase, latitude, human lactase phenotype, and IBD 



51. Ladas S, Papanikos J, Arapakis G. Lactose malabsorption in Greek 
adults: correlation of small bowel transit time with the severity of 
lactose intolerance. Gut. 1982;23(1 1):968-973. 

52. Sadre M, Karbasi K. Lactose intolerance in Iran. Am J Clin Nutr. 
1979;32(9):1948-1954. 

53. Nose O, Iida Y, Kai H, Harada T, Ogawa M, Yabuuchi H. Breath 
hydrogen test for detecting lactose malabsorption in infants and children. 
Prevalence of lactose malabsorption in Japanese children and adults. 
ArchDis Child. 1979;54(6):436-440. 

54. Nasrallah SM. Lactose intolerance in the Lebanese population and in 
Mediterranean lymphoma. Am J Clin Nutr. 1979;32(10):1994-1996. 

55. Aswami MZ, Seppo L, Vapaatalo H, Korpela R. Hypolactasia and 
lactose intolerance among three ethnic groups of Malaysia. Indian J 
Med Res. 2006;124(6):697-704. 

56. Senewiratne B, Thambipillai S, Pereira H. Intestinal lactase deficiency 
in Ceylon (Sri Lanka). Gastroenterology. 1977;72(6): 1257-1259. 

57. Filali A, Hassine B, Dhouib H, Matri S, BenAmmarA, Garoui H. Study 
of malabsorption of lactose by the hydrogen breath test in a population of 
70 Tunisian adults. Gastroenterol Clin Biol. 1987; 1 1(8-9):554-557. 

58. Flatz G, Henze HJ, Palabiyikoglu E, Dagalp K, TurkkanT. Distribution 
of the adult lactase phenotypes in Turkey. Trop Geogr Med. 1986;38(3): 
255-258. 

59. Bolin TD, Morrison RM, Steel JE, Davis AE. Lactose intolerance in 
Australia. MedJAust. 1970; 1(26): 1289-1292. 

60. Lember M, Torniainen S, Kull M, et al. Lactase non-persistence and 
milk consumption in Estonia. World J Gastroenterol. 2006; 12(45): 
7329-7331. 

61. Abbott WG, Tasman-Jones C. Incidence of acquired primary 
hypolactasia in three New Zealand racial groups. N Z Med J. 
1985;98(776):228-229. 

62. Farup PG, Monsbakken KW, Vandvik PO. Lactose malabsorption in a 
population with irritable bowel syndrome: prevalence and symptoms. 
A case-control study. Scand J Gastroenterol. 2004;39(7):645-649. 

63. Kuokkanen M, Butzow R, Rasinpera H, et al. Lactase and ovarian 
carcinoma risk in Finland, Poland and Sweden. Int J Cancer. 
2005;117(l):90-94. 

64. Mattar R, Monteiro MS, Villares CA, Santos AF, Silva JM, Carrilho FJ. 
Frequency of LCT-13910C.T single nucleotide polymorphism 
associated with adult-type hypolactasia/lactase persistence among 
Brazilians of different ethnic groups. Nutr J. 2009;8:46. 

65. Czeizel A, Flatz G, Flatz SD. Prevalence of primary adult 
lactose malabsorption in Hungary. Hum Genet. 1983;64(4): 
398^*01. 

66. Babu J, Kumar S, Babu P, Prasad JH, Ghoshal UC. Frequency of lactose 
malabsorption among healthy southern and northern Indian populations 
by genetic analysis and lactose hydrogen breath and tolerance tests. Am 
J Clin Nutr. 20 10;9 1 ( 1 ): 140-146. 

67. Al-Sanae H, Saldanha W, Sugathan TN, Majid Molla A. Comparison 
of lactose intolerance in healthy Kuwaiti and Asian volunteers. Med 
Princ Pract. 2003 ; 12(3): 1 60- 1 63 . 

68. Khabarova YA, Torniainen ST, Nurmi HA, Jarvela IE, Isokoski MK, 
Mattila KJ. Prevalence of lactase persistent/non-persistent genotypes 
and milk consumption in young population in north-west Russia. World 
J Gastroenterol. 2009; 15(1 5): 1 849-1 853. 

69. Socha J, Ksiazyk J, Flatz G, Flatz SD. Prevalence of primary 
adult lactose malabsorption in Poland. Ann Hum Biol. 1 984; 1 1 (4): 
311-316. 

70. Maggi R, Sayagues B, Fernandez A, et al. Lactose malabsorption and 
intolerance in Uruguayan population by breath hydrogen test (H2). 
JPediatr Gastroenterol Nutr. 1987;6(3):373-376. 

71. Ng SC, Bernstein CN, Vatn MH, et al. Geographical variability 
and environmental risk factors in inflammatory bowel disease. Gut. 
2013;62(4):630-649. 

72. Nerich V, Monnet E, Weill A, et al. Fine scale geographic variations 
of inflammatory bowel disease in France: correlation with socio- 
economic and house equipment variables. Inflamm Bowel Dis. 
2010;16(5):813-821. 



73. Armitage EL, Aldhous MC, Anderson N, et al. Incidence of juvenile- 
onset Crohn's disease in Scotland: association with northern latitude 
and affluence. Gastroenterology. 2004;127(4):1051-1057. 

74. Sonnenberg A, McCarty DJ, Jacobsen SJ. Geographic variation of 
inflammatory bowel disease within the United States. Gastroenterology. 
1991;100(1):143-149. 

75 . Khalili H, Huang ES, Ananthakrishnan AN, et al. Geographical variation 
and incidence of inflammatory bowel disease among US women. Gut. 
2012;61(12):1686-1692. 

76. Bernstein CN, Wajda A, Svenson LW, et al. The epidemiology of 
inflammatory bowel disease in Canada: a population based study. Am 
J Gastroenterol. 2006; 1 0 1(7): 1 559-1 568. 

77. Jussila A, Virta LJ, Salomaa V, Maki J, Jula A, Farkkila MA. High and 
increasing prevalence of inflammatory bowel disease in Finland with a 
clear north-south difference. J Crohns Colitis. 2013;7(7):e256-e262. 

78. Sonnenberg A. Time trends of mortality from Crohn's disease and 
ulcerative colitis. Int J Epidemiol. 2007;36(4):890-899. 

79. Garg M, Lubel JS, Sparrow MP, Holt SG, Gibson PR. Review 
article: vitamin D and inflammatory bowel disease - established 
concepts and future directions. Aliment Pharmacol Ther. 2012;36(4): 
324-344. 

80. Jorgensen SP, Agnholt J, Glerup H, et al. Clinical trial: vitamin D3 
treatment in Crohn's disease - a randomized double-blind placebo- 
controlled study. Aliment Pharmacol Ther. 2010;32(3):377-383. 

81. Cantorna MT, Mahon BD. Mounting evidence for vitamin D as an 
environmental factor affecting autoimmune disease prevalence. Exp 
Biol Med (Maywood). 2004;229(1 1):1 136-1 142. 

82. Holick MF, Chen TC, Lu Z, Sauter E. Vitamin D and skin physiology: 
a D-lightful story. J Bone Miner Res. 2007;22 Suppl 2:V28-V33. 

83. Aamodt G, Bengtson MB, Vatn MH. Can temperature explain the 
latitudinal gradient of ulcerative colitis? Cohort of Norway. BMC Public 
Health. 2013;13:530. 

84. Murdoch TB, Bernstein CN, El-Gabalawy H, et al. Prevalence of genetic 
variants associated with inflammatory bowel disease in a healthy First 
Nations cohort. CMAJ. 2012;184(8):E435-E441. 

85. Shapira M, Tamir A. Crohn's disease in the Kinneret sub-district, Israel, 
1960-1990. Incidence and prevalence in different ethnic subgroups. 
Eur J Epidemiol. 1 994; 1 0(2):23 1 -233 . 

86. Gearry RB, LeaRA, Roberts RL, Chambers GK, Barclay ML, Kennedy MA. 
CARD 15 allele frequency differences in New Zealand Maori: ancestry 
specific susceptibility to Crohn's disease in New Zealand? Gut. 
2006;55(4):580. 

87. Flatz G, Rotthauwe HW. Lactose nutrition and natural selection. Lancet. 
1973;2(7820):76-77. 

88. Heyer E, Brazier L, Segurel L, et al. Lactase persistence in central Asia: 
phenotype, genotype, and evolution. Hum Biol. 201 1 ;83(3):379— 392. 

89. Gallego Romero I, Basu Mallick C, Liebert A, et al. Herders of 
Indian and European cattle share their predominant allele for lactase 
persistence. Mol Biol Evol. 2012;29(l):249-260. 

90. Szilagyi A, Leighton H, Burstein B, Shrier I. Significant positive 
correlation between sunshine and lactase nonpersistence in Europe 
may implicate both in similarly altering risks for some diseases. Nutr 
Cancer. 201 1;63(7):991-999. 

91. Szilagyi A. Lactose: a potential prebiotic. Aliment Pharmacol Ther. 
2002;16(9):1591-1602. 

92. Szilagyi A, Nathwani U, Vfnokuroff C, Correa JA, Shrier I. The effect 
of lactose maldigestion on the relationship between dairy food intake 
and colorectal cancer: a systematic review. Nutr Cancer. 2006;55(2): 
141-150. 

93. Spooren CE, Pierik MJ, Zeegers MP, Feskens EJ, Masclee AA, 
Jonkers DM. Review article: the association of diet with onset and 
relapse in patients with inflammatory bowel disease. Aliment Pharmacol 
Ther. 2013;38(10):1 172-1 187. 

94. Abubakar I, Myhill DJ, Hart AR, et al. A case control study of drinking 
water and dairy products in Crohn's disease: further investigation of 
the possible role of Mycobacterium avium paratuberculosis. Am J 
Epidemiol. 2007;165(7):776-783. 



Clinical Epidemiology 2014:6 



submit your manuscript | www.dovepress.com 
Dovepress 



197 



Szilagyi et al 



Dovepress 



95. Alfonso-Sanchez MA, Perez-Miranda AM, Garcia-Obregon S, Pena JA. 
An evolutionary approach to the high frequency of the delta 
F508 CFTR mutation in European populations. Med Hypotheses. 
2010;74(6):989-992. 

96. Gasche C, Nemeth M, Grundtner P, Willheim-Polli C, Ferenci P, 
Schwarzenbacher R. Evolution of Crohn's disease-associated Nod2 
mutations. Immunogenetics . 2008;60(2):1 15-120. 

97. Wang TT, Dabbas B, Laperriere D, et al. Direct and indirect induction 
by 1,25-dihydroxyvitamin D3 of the NOD2/CARD15-defensin 02 
innate immune pathway defective in Crohn's disease. J Biol Chem. 
2010;285(4):2227-2231. 

98. Modiano G, Ciminellui BM, Pignatti PF. Cystic fibrosis and lactase 
persistence: a possible correlation. Eur J Hum Genet. 2007;15(3): 
255-259. 

99. Bresso F, Askling J, Astegiano M, et al. Potential role for the common 
cystic fibrosis AF508 mutation in Crohn's disease. Inflamm Bowel Dis. 
2007;13(5):531-536. 

1 00. Bahmanyar S, Ekbom A, Askling J, Johanneson M, Montgomery SM. 
Cystic fibrosis gene mutations and gastrointestinal disease. J Cyst 
Fibros. 2010;9(4):288-291. 

101. Buning C, Ockenga JK, Jurga J, et al. The C/C(-13910) and G/ 
G9-22018) genotypes for adult-type hypolactasia are not associated 
with inflammatory bowel disease. Scand J Gastroenterol. 2003;38(5): 
538-542. 



102. Elguezabal N, Chamorro S, Molina E, et al. Lactase persistence, 
NOD2 status and Mycobacterium avium subsp. Paratuberculosis 
infection associations to inflammatory bowel disease. Gut Pathog. 
2012;4(1):6. 

103. Nolan DJ, Han DY, Lam WJ, et al. Genetic adult lactase persistence is 
associated with risk of Crohn's disease in a New Zealand population. 
BMC Res Notes. 2010;3:339. 

104. Juste RA, Elguezabal N, Chamorro S, et al. Associations between 
CARD 1 5 polymorphisms, MAPDNA in blood and lactase persistence 
in a Crohn's disease case-control study in north Spain. Proceedings 
of the Tenth International Colloquium on Paratuberculosis. 2009; 
269(2):9-14. 

1 05. Lees CW, Barrett JC, Parkes M, Satsangi J. New IBD genetics: common 
pathways with other diseases. Gut. 201 1;60(12): 1739-1753. 

106. Bernstein CN, Kraut A, Blanchard JF, Rawsthorne P, Yu N, Walld R. The 
relationship between inflammatory bowel disease and socioeconomic 
variables. Am J Gastroenterol. 2001;96(7):21 17-2125. 

107. Ross AG, Olds GR, Cripps AW, Farrar JJ, McManus DP. 
Enteropathogens and chronic illness in returning travelers. N Engl J 
Med. 2013;368(19):1817-1825. 

108. Okada H, Kuhn C, Feillet H, Bach J-F. The 'hygiene hypothesis' for 
autoimmune and allergic diseases: an update. Clin Exp Immunol. 
2010;160(l):l-9. 



Dovepress 



Clinical Epidemiology 

Publish your work in this journal 

Clinical Epidemiology is an international, peer-reviewed, open access 
journal focusing on disease and drug epidemiology, identification of 
risk factors and screening procedures to develop optimal preventative 
initiatives and programs. Specific topics include: diagnosis, prognosis, 
treatment, screening, prevention, risk factor modification, systematic 

Submit your manuscript here: http://www.dovepress.com/clinical-epidemiology-journal 



reviews, risk & safety of medical interventions, epidemiology & bio- 
statical methods, evaluation of guidelines, translational medicine, health 
policies & economic evaluations. The manuscript management system 
is completely online and includes a very quick and fair peer-review 
system, which is all easy to use. 



submit your manuscript | www.dovepress.com 

Dovepress 



Clinical Epidemiology 2014:6 



